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Finite -element methods along with Monte 
Carlo simulations were used to design a 
magnetic storage device for ultracold 
neutrons (UCN) to measure their lifetime. 
A setup was determined which should 
make it possible to confine UCN with 
negligible losses and detect the protons 
emerging from p-decay with high 
efiiciency: stacked superconducting 
solenoids create the magnetic storage field, 
an electrostatic extraction field inside 
the storage volume assures high proton 
collection efficiency. Alongside with 
the optimization of the magnetic and 
electrostatic design, the properties of the 
trap were investigated through extensive 
Monte Carlo simulation. 
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1. Introduction 

The neutron lifetime \ plays a vital role in under- 
standing big bang cosmology: it influences the relative 
abundance of primordial helium in the early universe. Tn 
also opens the way to determine the coupling constants 
of the weak interaction and hence the element V^^ of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix precisely 
Latest experimental results indicate this matrix to 
deviate from unitarity by about 3 standard deviations 
and therefore challenge the three-generation Standard 
Model [1]. 

The most precise measurements of T^ were per- 
formed by storing UCN in material bottles. There are, 
however, significant losses during the many wall 
collisions whose nature is not yet fully understood. 
Therefore systematical errors cannot be decreased 
much below their present values. Magnetic storage has 
recently been proven to be a viable alternative [2]. This 
publication presents a new experiment, designed for 
new generation UCN sources as the one planned for the 



research reactor FRM II at Garching [3]. It utilizes 
magnetic storage in combination with real-time proton 
detection. Furthermore, the results of Monte Carlo 
calculations for the proposed setup are shown. 

2. Description of tlie Experiment 

Magnetic storage is based on the force F on a magnet- 
ic moment fi that is exerted in an inhomogeneous mag- 
netic field with flux density B 



F = V{ixB) 



(1) 



Only particles with one orientation of the spin towards 
B can be stored, hence reorientation (also called 
depolarization) has to be avoided to assure loss-less 
neutron storage. To reach this goal, the adiabatic 
condition has to be fullfilled: changes B of the magnetic 
field (especially rotations of the field vector) seen by 
the moving neutron normalized to the absolute value of 
the magnetic flux density B have to be much slower 



357 



Volume 110, Number 4, July-August 2005 

Journal of Research of the National Institute of Standards and Technology 



than the neutron Larmor precession in the magnetic 
field o^L 



B 



«0), 



(2) 



For the proposed setup, the volume between two 
nested cylinders made from magnetic multipole fields 
is used to store the UCN (Fig. 1). This storage field is 
produced by stacked superconducting coils of alternat- 
ing current direction with the gravitational field playing 
the role of the upper lid of the bottle. The dimensions 
of 1 .2 m height and 0.5 m outer radius result in a storage 




Fig. 1. Cutaway view of the experimental arrangement. See text for 
further explanations. 



volume of V ~ 0.8 m"^. In the center of the setup an 
additional current carrying rod creates an azimuthal 
magnetic field, which is always perpendicular to the 
storage field and therefore helps keep the flux density 
in the whole storage volume above a critical value, 
below which reorientation of the spin may occur 
[see Eq. (2)]. Filling and emptying of UCN is realized 
through a slit in the outer bottom comer while the 
current in the storage coils is low enough to let neutrons 
enter the trap. 

Tn shall be measured by real-time detection of decay 
protons as well as by counting the integral number of 
neutrons using different storage times. The proton 
detector, a Csl scintillator, is situated on top of the 
storage volume; the decay protons are accelerated and 
focused onto it through a potential difference in the 
storage volume and an additional focusing coil around 
the detector. Depolarized neutrons shall also be moni- 
tored; they are not stored magnetically, but can still be 
collected when the inner walls of the trap are covered 
with a neutron reflecting material as, e.g., diamond-like 
carbon. Using a new high-density UCN source at the 
FRM II [3] it is possible to store up to 10^ neutrons 
per cycle. Hence, the measuring time to get sufficient 
statistics is short and many runs with different condi- 
tions can be realized to investigate possible systematic 
errors. We thus envisage a relative accuracy for 
T.oflO^. 

3. Monte Carlo Simulation 

The behavior of all particles involved in neutron 
decay in the experiment is dominated by the influence 
of the magnetic field in the case of neutrons or com- 
bined magnetic and electrostatic fields for the charged 
decay particles (neutrinos may be neglected). The fields 
(see Fig. 2) were calculated using three different finite 
element method programs [4-6], which produced the 
same results. These were then used as an input for 
Monte Carlo simulations. A Runge-Kutta algorithm 
with adaptive step-size control was used to calculate 
the trajectories of protons and neutrons [7]. Example 
trajectories for a stored neutron, a depolarized neutron 
and several protons are displayed in Fig. 3 . 

It was shown that neutrons of less than 120 neV kinet- 
ic energy can be stored loss-less inside the trap. A filling 
and emptying time for neutrons of less than 50 s could be 
confirmed as sufficient. Depolarized neutrons, arising 
from spin flip of the stored UCN, but also present after 
filling, are removed after about 25 s. When installing a 
neutron detector below the storage volume, they can be 
detected with an efficiency of around 60 %. 
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Fig. 2. Contour plot of the magnetic flux density |5|, including representative lines of constant \B\. 



Another critical issue that could change the meas- 
ured lifetime was also investigated: neutrons of kinetic 
energy E]^^ higher than 120 neV have a shorter storage 
time, most of them (96 %) are gone after 100 s, but the 
ones close to 120 neV have to be cleaned by inserting 
an absorbing material into the storage volume at a 
position only accessible to UCN of ^ki^ > 120 neV. 

Using Vladimirsky's approach [8], the probability 
for UCN depolarization was found to depend on the 
current in the center rod: for a total current of 3 kA this 
probability (3 x 10"^) is already well below that required 
to reach the desired lifetime accuracy 



Protons emerging from neutron decay in the storage 
volume can be collected at the detector with an 
efficiency close to 80 %, their kinetic energies there are 
in the range from 30 keV to 40 keV and hence high 
enough for detection. The background of electrons 
at the scintillator was determined to be manageable, 
as only 2 % of all decay electrons arrive at its posi- 
tion. Furthermore, less than 10"^ of them deposit 
enough energy in the scintillator to be confiised 
with the proton signal. This is again well below our 
accuracy goal. 
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Fig. 3. r-z components of example trajectories: (left) a stored neutron (r<25s) — (center) a depolarized neutron (r<10s) — (right) five 
decay protons. 



4. Summary and Status 

The neutron lifetime is a fundamental and important 
constant of nature. Using magnetic storage of UCN 
it may be measured with high precision. The feasibility 
of the proposed design from the physics point of view 
has been confirmed through Monte Carlo simulation of 
the involved particles. The realization is on the way the 
setup is being optimized for technical feasibility and first 
assemblies could start later this year. 
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